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Abstract 

First correction to the high-energy asymptotics of the total e'^e" photoproduction 
cross section in the electric field of a heavy atom is obtained with the exact account 
of this field. The consideration is based on the use of the quasiclassical electron Green 
function in an external electric field. The influence of screening on the Coulomb 
corrections is examined in the leading approximation. It turns out that the high- 
energy asymptotics of the corresponding correction is independent of the photon 
energy. The detailed comparison of our results with experimental data is performed. 
This comparison has justified the analytical result and allowed us to elaborate a 
simple ansatz for the next-to-leading correction. Using this ansatz, good agreement 
with the experimental data is obtained for photon energies above a few MeV . In the 
region where both produced particles are relativistic, the corrections to the high- 
energy asymptotics of the electron (positron) spectrum are obtained. In addition, 
analogous corrections to the bremsstrahlung spectrum are derived starting from the 
corresponding results for pair production. 

Key words: e^e^ photoproduction, bremsstrahlung. Coulomb corrections, 
screening 
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I Introduction 



Knowledge of the phot oabsorpt ion cross sections is very important in various 
applications, see, e.g., ( Hubbelj 2(]00). The relevant processes are the atomic 



Email addresses: R.N.Lee@inp.nsk.su (R.N. Lee), A.I.Milstein@inp.nsk.su 
(A. I. Milstein), V.M.Strakhovenko@inp.nsk.su (V.M. Strakhovenko). 



Preprint submitted to Elsevier Science 



7 February 2008 



photoeffect, nuclear photoabsorption, incoherent and coherent photon scatter- 
ing and e~^e~ pair production. In the coherent processes, by definition, there 
is no excitation or ionization of an atom. The high-accuracy estimation of the 
corresponding cross sections is required. They have different dependence on 
the photon energy uj. At u; > lOMeV, the cro ss section of e^e^ pair pro- 



duction becomes dominant ( Hubbell et al. . IQSOi ). The coherent contribution 



o'coh to the pair production cross section is roughly Z times larger than the 
incoherent one {Z is the atomic number), thereby being the most important 
for heavy atoms. Just the coherent pair production is considered below. 



The theoretical and exp erimental investigati on of the coherent pair production 



has a long history, see ( Hubbell et al. . 1980l ). In the Born approximation, the 



cross section ctr is known for arbitrary photon energy (jBethe and Heitlei , 
1934t iRacahL Il934^ . The account of the effect of screening is straightforward 
in this approximat i on an d can be easily performed if the atomic form factor is 
known ( Jost et al. . IQSOl l. For heavy atoms it is necessary to take into account 
the Coulomb corrections ac, 



(^coh = CTb + CTC ■ 



These corrections are higher order terms of the perturbation theory with re- 
spect to the atomic field. The magnitude of ac depends on u and the param- 
eter Za {a = 1/137 is the fine-structur e constant). T h e form al expression for 
ac, exact in Za and u, was derived by 0verb0 et al. ( 1968[ ). This expression 



has a very complicated form causing severe difficulties in com putations. The 



diffic ulties grow as u increases, so that numerical results in (|0verb(3 et al, 
19681) were obtained only for cu < 5MeV. 



For the high-energy region uj ^ m [m is the electron mass), the considera- 
tion is greatly simplified. As a result, a rather simple form was obtained in 
( Bethe and Maximonl . 119541 iDavies et al.l . Il954^ for the Coulomb corrections 
in the leading approximation with respect to m/u. However, the theoreti- 
cal description of the Coulomb corrections at intermediate photon energies 
(5 -ir lOOMeV^) has not been completed. At present, all es timates of err; in this 
region are based on the "bridging" expression derived by 0verbd (1977). This 
expression is actually an extrapolation of the results obtained for uj < 5MeV. 
It is based on some assumptions on the form of the asymptotic expansion of 
ac at high photon energy. It is commonly believed that the "bridging" ex- 
pression has an accuracy providing the maximum error in acoh of the order of 
a few tens of percent. 



Here we develop a description of e~^e -pair production at intermediate photon 
energies by deriving the next-to-leading term of the high-energy expansion of 
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(7c- First we consider a pure Coulomb field and represent ac in the form 



The term o"^'' has the form (m/co')"5'*^"''(lnti;/m), where 5'*^"''(a;) is some poly- 
nomial. The a;-independent term cr^'' corresponds to the result ofE avies et alJ 



( 1954 ). In the present paper we derive the term a^^. It turns out tha t .S**-^-* is 



a;-inde pendent in contrast to a second-degree polynomial suggested bv l0ver 
( 19771 ). We propose a new ansatz for , which provides a good agreement 



with available experimental data for uj > 5MeV. 

The high-energy expansion of the Coulomb corrections to the spectrum has the 
same form as (2). In the region e± S> m, we derive the term da^^ /dx, where 
E- and are the electron and positron energy, respectively, x = E-joj. The 
term da^^ /dx may turn important, e.g., for description of the development of 
electromagnetic showers in a medium. The correction found is antisymmetric 
with respect to the permutation ^ e_ and does not contribute to the total 
cross section. In fact, originates from two energy regions ~ m and 
E- ~ m, where the spectrum is not known. However, we emph asize that it 



differs drastically from the result obtained bv lDavies et al.l ((195J) for e-t ^ m 



if the latter is formally applied at £_ ~ m or ~ m. 

For the first time, we estimate the effect of screening on ac at u; ^ m. In the 
leading approximation, we find the corresponding correction (yc"'\ which is 
cu-independent similar to cx^'*. So, for the atomic field, cr^'^'^'' should be added 
to the right-hand side of Eq. (2). The screening correction to the spectrum is 
also obtained. 

Tedious and cumbersome calculations have been performed to derive our for- 
mulas. All technical details are omitted here and will be presented elsewhere. 



II General discussion 



The cross section of e^e pair production by a photon in an external field 
reads 

d(Tcoh= ,^ dpdq5{uj - E+- eJ)\M\^ , (3) 



where E- — Ep — \pp^^^\^'rr? ^ e+ = Sq, and p, q are the electron and positron 
momenta, respectively. The matrix element M has the form 

M = j drij^^\r)e'ilj^~\r)exp{ikr) . (4) 



Here and tpl~^ are positive-energy and negative-energy solutions of the 
Dirac equation in the external field, is the photon polarization 4- vector, k is 
the photon momentum. It is convenient to study various processes in external 
fields using the Green function G(r2, rije) of the Dirac equation in this field. 
This Green function can be represented in the form 



E — En +iO 



+ 



J 



dp 

(2^ 



E — Ep +iO 



+ 



E + Ep — iO 



(5) 



where ipii'^ is the discrete-spectrum wave function, e„ is the corresponding 
binding energy. The regularization of denominators in (5) corresponds to the 
Feynman rule. From (5), 



I dn^^^^\rM^^\r,)^z^^5G{r,,r,\Ep), (6) 

where flp is the solid angle of p, and SG — G — G. The function G is obtained 
from (5) by the replacement iO ^ —iO. 

Taking the integrals over Jlp and flq in (3), we obtain the electron spectrum, 
which is the cross section differential with respect to the electron energy £_ . 
Using the relations (6), we express this spectrum via the Green functions: 



-^ = ^11 dr.dr^e'^^Sp {5G{r,,r2\E.) e6G{r2,r,\ - e+) t} , (7) 
where e+ — u — e^ and r — r2 — ri. 

Due to the optical theorem, the process of pair production is related to the 
process of Delbriick scattering (coherent scattering of a photon in the electric 
field of an atom via virtual electron-positron pairs). At zero scattering angle, 
the amphtude of Delbriick scattering reads 
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MD = 2ia jde jj dri dr2e''^'^ Sp {G{ri,r2\e) eG{r2,ri\e - u) e*} 



It follows from Eqs. (7), (8) and the analytical properties of the Green function 
that 

-Im. Mr, = acoh + (^bf ■ (9) 



Here 



(Jbf = f [ dridr2e'^'^Yl Sp {pn{ri,r2) e6G{r2,ri\en - uj) e*} , 

LO J J „ 

Pn(n, r2) = lim {e - £:„)G(ri, rsje). (10) 



The quantity abj coincides with the total cross section of the so-called bound- 
free pair production when an electron is produced in a bound state. In fact, due 
to the Pauli principle, there is no bound-free pair production on neutral atoms. 
Nevertheless, the term a^f should be kept in the r.h.s. of (9). In a Coulomb 
field, the total cross section af,/ was obtained in ( Milstein and Strakhovenkol . 



Il993r ) for 3> m. In this limit, abj oc 1/muj and should be taken into ac- 
count when using the relation (9) for the calculation of the corrections to acoh 
from. The main contribution to (T hf comes from the low-lying bound states 
(iMilstein and Strakhovenkl Il99.^ when screening can be negle cted. So, in 



(9) we can use abj obtained in ( Milstein and Strakhovenkol . 1993[ l. 



Ill Coulomb corrections to the spectrum 



In this section we consider the Coulomb corrections to the spectrum, dac/dx, 
for £± ^ m taking into account terms of the order m/e±. Within this accuracy, 
the spectrum is determined by the region of small angles between vectors 
p, g, and k in (3). In this case the angular momenta of both particles are 
large, providing the applicability of the quasiclassical approximation. Thus we 
can use the quasiclassical G reen function in (7). For a pure Coulomb field, 
this function was found in ( Milstein and Strakhovenkol . Il983h : for arbitrary 



localized field, it was obtained in (jhee et al.l . 20001 ). In the latter paper, the 
first correction to the quasiclassical Green function was derived as well. 



According to lDavies et al.l (|195J), the higher order terms of the perturbation 
theory with respect to the external field (Coulomb corrections) are not se- 
riously modified by screening. However, this question has not been studied 
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quantitatively so far. The influence of screening on Coulomb corrections is in- 
vestigated in detail in Section V. In the present Section we calculate dac/de_ 
in a pure Coulomb field. 



Using the results of iLee et al. (|2nnnl) . we obtain from (7) 



1-^xil-x)) fiZa) 



^ ^ ' 1 - -x) Re^(Za) 



8x(l-x)cu V 2 

. , . ^ , . T(l-iZa)T(l/2 + iZa) 

x = e-/uj , (To = a(2'a)^/m^ , (11) 

where il){t) = dhir{t)/dt, C = 0.577... is the Euler constant. In (11), the term 
PC f{Za) corresponds to the leading approximation da^^ /dx (jPavies et al 



1954 ) ■ the term oc Reg{Za) is the first correction da^^^'/dx. In contrast to 



the leading term, this correction is antisymmetric with respect to the permu- 
tation ^ £_ (or X ^ 1 — x) and, therefore, does not contribute to the 
total cross section. Besides, the correction is an odd function of Za due to 
the charge-parity conservation and the antisymmetry mentioned above. The 
antisymmetric contribution enhances the production of electrons at x < 1/2 
and suppresses it at x > 1/2. Evidently, the opposite situation occurs for 
positrons. Qualitatively, such a behavior of the spec trum takes place for any 
u being the most pronounced at low photon energy ( 0verb(Z) et al. . 1968| ). At 



intermediate photon energies, the spectrum (11) essentially differs from that 
given by the leading approximation. We illustrate this statement in Fig. 1, 
where a^^dac I dx with correction (solid line) and without correction (dashed 
hne) are plotted for Z = 82 and = 50 MeV. 

Due to the antisymmetry of da^^ /dx at e± 3> m, the term a^^ in the total 
cross section may originate only from the energy regions S- ^ m and e+ = 
uj — E- m. The quasiclassical approximation can not be used directly in 
these regions, and another approach is needed to calculate the spectrum. We 
are going to do this elsewhere. However, for the total cross section, it is possible 
to overcome this difficulty by means of dispersion relations (see Section IV). 



As known (see, e.g., Berestetski et al. . 1982| ). the spectrum of bremsstrahlung 



can be obtained from the spectrum of pair production. This can be performed 
by means of the substitution e+ — * — e , ^ — c<j', and dx ydy, where 
y = uj'/e, ijj' is the energy of an emitted photon, e is the initial electron en- 
ergy. Using (11), we obtain for the Coulomb corrections to the bremsstrahlung 
spectrum 
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Fig. 1. The dependence of aQ dac/dx on x, see (11), for Z = 82, to = bOMeV. 
Dashed curve: leading approximation; solid curve: first correction is taken into ac- 
count. 



' dy 



-4ao 

7r'^(2 — y)m 



y) fiZa) 



y^ + -,ii 



y)] Reg{Za] 



(12) 



This formula describes bremsstrahlung from electrons. For positrons, it is nec- 
essary to ch ange the sign of Za in (1 2). Our result (12) coincides with that 
obtained in ( Baier and Katkov . 1976t ) if the obvious mistake in the latter is 
corrected by changing 



1 
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m m 
— + 

e e 



UJ' 



(2 - y)m 
2(l-2/)£ 



in Eq.(22) of ()Ra,ier and Katkovl . Il 97(tI ) 



IV Coulomb corrections to the total cross section 



For ^ m, we derive cx^^ using the relation (9). In the leading approximation, 
the Coulomb correction s to the cross section of pair production were obtained 



in (jPavies et al.L Il954l ). Using this result and dispersion relations, the corre- 
sponding correcti on, Mnr^, to the forward Delbriick scattering amplitude Md 
were obtained in (jRohrlichL Il957l ). They read: 
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a. 



(0) 

c 



(13) 



where cxo and f{Za) are defined in (11). 

Using the results of Ihee et aL ( 200(1 ) for the quasiclassical Green functions, 
we find within logarithmic accuracy 



m m 



(14) 



the function g{Za) is defined in (11). The logarithm appears as a result of 
integration in (8) over e in the region where e ^ m and u — e ^ m, thereby 
the quasiclassical approximation is applicable. Unlike lie Mj^q, the quantity 
ImM^j^ is determined by the regions of integration over e, where e ~ m 
and uj — 6 ^ m, and, therefore, the quasiclassical approximation is invalid. 
Nevertheless, this quantity, whi ch is related to a^^ (9), can be obtained from 
the dispersion relation for Md ()Rohrlichl . 119571 ) 



OO 

ReM„M = -.^p/-;^;^-i^. (15) 







Using this relation, it can be easily checked that the high-energy asymptotics 
(14) unambiguously corresponds to the a;-independent high-energy asymp- 
totics 

ImMa = - "^^"^'"'^"^^(^"\ (16) 

2m 



Subst ituting (16) into (9) and using abf from ( Milstein and Strakhovenkol . 
I993I ) in the form 



a,f= AiraoiZaf h{Za)- , (17) 

UJ 

we have for ctq^ 



(1) 



IT 



— lmg{Za) + AiT{Zaf fi{Za) 



m 



(1^ 



The function fi{Za) is plotted in Fig. 2. 
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Fig. 2. The quantity /i as a function of Z 

The quantity {uj j vnb)o^^ j g^q is shown in Fig. 3 (sohd curve). It is seen that this 
ratio is numerically large for any Z . Therefore, the term a^'' gives a significant 
contribution to gq for intermediate photon energies. Dashed curve in Fig. 3 
gives the same ratio when ct;,/ in (18) is omitted. It is seen that the relative 
contribution of the term oc fi{Za) in (18) is numerically small. 
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Z 

Fig. 3. The quantity (a;/m)o"^''/cj^^ as a function of Z (solid curve). Dashed curve 
corresponds to the same quantity without the contribution of the bound-free pair 
production. 
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V Screening corrections 



In two previous Sections the cross section of e^e~ pair production has been 
considered for a pure Coulomb field. The difference, 5V{r), between an atomic 
potential and a Coulomb potential of a nucleus leads to the modification of this 
cross section known as the effect of scr eening. In the Bo rn approximation, this 
effect was studied long ago (see, e.g., Ijost et al.l . ll95nh . Let us consider now 



o"^'^'^'' characterizing the influence of screening on the Coulomb corrections. 
Recollect that the Coulomb corrections denote the higher-order terms of the 
perturbation theory with respect to th e atomic field. So far it was only known 



that the correction a^^'^^'' is not large (|Davies et al.L Il954( ). Here we consider 



this issue quantitatively. The Coulomb corrections are determined by distances 
r ~ 1/m where the difference 5V{r) is small. In our calculation of (Jc'^'^\ we 
retain the linear term of expansion in powers of 5V{r). We represent 5V{r) as 



5V{r) 



dQ 

(2^ 



FiQ) 



AirZa 



(19) 
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Fig. 4. The ratio a^^^^ /(Tq^ as a function of Z 

where F(Q) is the atornic electron form factor. Using the results obtaine d in 
fjLee and Milsteinl . Il99,4 iKatkov and Strakhovenkol . boOlULee et aLlboOOl ) for 

arbitrary atomic potential, we have for the correction to the spectrum 



da. 



(scr) 
C 



32 



dx 



-aom 



dQ 
Q 



dr 



sinh r 



sin(2ZQ;r) 
2Z^ 



10 



2-K 



dip r ^ 



J 2tt 



a) 



e^R{fi+, a) — e '^R{fi-, a 
1 



-18-a(/i-3)| 

8m^ e^"^ sinh^ r 

/i± = 1 + 



(20) 



18 - 6/i + a(/i^ + 2/i - 3) 



In 



(5^(coshr + cos(p) 
Integrating over x, we obtain 



a = 6x(l — x) . 



(21) 



a, 



(scr) 32 _ ^2 



C 



dQ 



dr 



sinh r 



sin(2Zar) 
2Z^^ 



2lT 



dip r ^ 



27r 



e^i?(/i+, 1) -e"^i?(/i_, 1) 



(22) 



Similar to a'^\ this correction is tu-independent. Shown in Fig. 4 is the Z- 
dependence of the ratio an^^^ / qff | calcul ated with the use of the form factors 
taken from ( Hubbell and 0verbeil . 1979| ). As seen from Fig. 4, this ratio is 
approximately described by the linear function —5.4 • 10~^ ■ Z. 

The corresponding correction to the bremsstrahlung spectrum is obtained from 
(20) by means of the same substitutions as in Section III. So that the quantity 
y~^da]}^^^'^ /dy is given by the right-hand side of (20) if we set a = 6{y — l)/?/^- 



(2) 



VI Comparison with experimental data and estimation of 



The most detailed and accurate experimental data have been obtained just in 
the region of intermediate photon energies. In this region, the first correction 
0"^'*, obtained above, becomes large, see Fig. 3, and the next term in the 
expansion (2) may be significant. Since it has not been calculated, we use for 
0"^ the following ansatz: 



(J. 



(2) 



Q — aQ[b\n{iu/2Tn) + c] [ — j 



(23) 



where b and c are some functi ons of Za. Th i s form follows from the arguments 
similar to those presented bv lDavies et al.l ()l954r ). 
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Shown in Fig. 5 is the quantity 



^ — — CTq K'^coh 

m 



(J. 



(0) 

c 



0", 



(scr) 
C 



(24) 



wher e the experimental cross section acoh for Bi is taken from ([Sherman et al 
1980l ). aB is the Born cross section calculated with screening taken into ac- 



(0) 

c 



and a^,'^^^ are given by (13) and (22), respectively. This quantity is 



count, a, 

fitted by the formula a + (m/to')[61n(ti;/2m) + c] (dashed curve) which corre- 
sponds to the sum {uj/m)aQ^{o'c^ + ctc^). The fitting parameters obtained by 
the linear regression method are a = 25.79, b = — 98.92, c = 2.43. The quantity 
a determines the asymptotics of the fit at u; ^ oo and should be compared 
with our result ath = /'^)<^o^ ctc"* — 25.66 (sohd line). So the value of ath is 
in a perfect agreement with that extracted from the experimental data. From 
Fig. 5, we conclude that the term a^'' gives a noticeable contribution to S. 
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uj {MeV) 

Fig. 5. The values of S (24) extracted from the experimental data for Bi together 
with the fit a+(m/u;)[61n(u;/2m)+c] (dashed cmve) and the theoretical asymptotics 
ath = ('^/"^)'^o'^'^c^ (solid line). Coefficients a, 6, and c are given in the text. 



Note that the coefficient c is small as compared to h. The sam e situation (a «: 
ath. \ b\ S> |c| ) takes place at fitting of experimental data for Ta ([Sherman et al. 
1980). So, we set c = in Eq. (23) and use h as the only fitting parameter 



Fitting the data for Bi with such ansatz, we obtain h = —96.63. A new curve 
practically coincides with that shown in Fig. 5. 



The values of S extracted from the experimental data for Ta dSherman et al. 



E3 and Pb (jRosenbhim et a1.l . 119521 : iGimm ajid Hubbe1ll . h978t ) are sh"^ 
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in Fig. 6 and Fig. 7, respectively. Dashed curves represent the results of the 
fit, solid lines correspond to the asymptotics = ath- 




5 10 15 20 25 30 

uj (MeV) 

Fig. 6. The values of S (24) extracted from the experimental data for Ta together 
with the fit ath + (m/a;)61n(u;/2m) (dashed curve), ath = 21.88, b = —82.82. The 
solid line represents the asymptotics S = ath- 
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UJ (MeV) 

Fi g. 7. The values of T, (24) extrac ted from the ex perimental data for Pb obtained 
bv * E,osenblum et al.l (squares) and lSherman et al. l (stars). The dashed curve corre- 
sponds to the fit ath + {rn/Lj)bln{io/2m), ath = 25.29, b = —95.95. The solid line 
represents the asymptotics S = ath- 

Within the accuracy of the fitting procedure, we obtain b = —3.78 ■ ath for all 
three elements considered. This fact hints at the same Z-dependence of the 
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coefficients ath and b, at least for heavy atoms. 



In the paper (jPavies et aLl |l95J) it was mentioned that all available ex- 
perimental data on Pb above 5MeV are well represented by the formula 
o'coh = (Tb + o"c + 11-8 • {m/uj)(ro. Poor knowledge of the nuclear photoabsorp- 
tion cross section at that time can not explain the large difference between the 
coefficient 11.8 and our result ath = 25.29. The main source of this difference 
is the neglecting of the term . It is evident from Fig. 7 that an attempt to 
fit E in the region below 18 MeV by the constant only, leads immediately to 
the appreciably smaller coefficient a than our value ath- 



It is interesting to compare our predicti ons for the Co ulomb corrections to 



of4°\ 



0verb(zJ (|l977t) . Shown in Figs. 8,9 



the total cross section with the results of 

is the ratio S = (acoh ~ o'B)/crr\ which is the Coulomb corrections in units 



(13). Our results are represented by solid curves, those of 10 verba are 
shown as dashed curves. The values of 5* extracted from the experimental 
data are also shown. The results for Bi are plotted in Fig. 8 with the ex- 



perimental data taken from ( Sherman et al. . Il980l ). The results for Pb are 



plott e d in Fig. 9 with the expe rimental data taken from ([Rosenblum et al 
19521: iGimm and Hubbell [l978l ). It is seen that the difference between our 
results and those of 0verbd is small at relatively low energies and becomes 



noticeable as u increases. According to our results, this difference tends to a 
constant cr^'^^V'^c'' ^.t ui ^ oo. As a whole, the experimental data are in a 



better agreement with our results than with those of l0ver 
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Fig. 9. Same as Fig. 8 but for Pb; experimental data from (jRosenblum et al 
Gimm and Hubbell llQT^ 'l. 
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VII Conclusion 



For the e^e~ photoproduction, we have calculated the leading correction (11) 
to the electron spectrum in the region e± ^ m. This contribution noticeably 
modifies the spectrum at intermediate photon energy. It turns out that the 
correction is antisymmetric with respect to the permutation e+ ^ and 
hence does not contribute to the total cross section. The leading correction to 
the total cross section, (Jq\ originates from two regions, ~ m and ~ m. 
We have obtained a^^ (18) using dispersion relations. In contrast to the form 
of the fit suggested bv l0verb0l (1977), the quantity does not contain any 
powers of ln(ti;/m). We have also performed the quantitative investigation of 
the influence of screening on the Coulomb corrections (20), (22). It is important 
that o"^'^'^'' does not vanish in the high-energy limit. We have suggested a form 
for the next-to- leading correction, cr^ , to the total cross section. Altogether, 
the corrections found allow one to represent well the available experimental 
data. 



Starting with the results obtained for the e^e~ photoproduction spectrum, we 
have obtained the corresponding correction to the bremsstrahlung spectrum 
as well. 
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